Subgap states in dirty superconductors and their effect on dephasing in Josephson 

qubits 
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We present a theory of the subgap tails of the density of states in a diffusive superconductor 
containing magnetic impurities. We show that the subgap tails have two contributions: one arising 
from mesoscopic gap fluctuations, previously discussed by Lamacraft and Simons, and the other 
associated to the long-wave fluctuations of the concentration of magnetic impurities. We study 
the latter both in small superconducting grains and in bulk systems [d = f , 2, 3], and establish the 
dimensionless parameter that controls which of the two contributions dominates the subgap tails. 
We observe that these contributions are related to each other by dimensional reduction. We apply 
the theory to estimate the effects of a weak concentration of magnetic impurities [« lp.p.m] on 
the phase coherence of Josephson qubits. We find that at these typical concentrations, magnetic 
impurities are relevant for the dephasing in large qubits, designed around a f fim scale, where they 
limit the quality factor to be Q < 10 4 — f0 5 . 

PACS numbers: 74.40.+k; 74.81.-g; 85.25.-j. 



I. INTRODUCTION. 

Recent experimental studies of disordered conductors 
demonstrated that even a weak concentration of mag- 
netic impurities may lead to important effects, especially 
in the context of dephasing and energy relaxation. In par- 
ticular, the experiments in mesoscopic Au and Cu wires 
have shown that the frequently observed saturation of 
the dephasing time t v at low temperatures^^ is likely 
due to a small unavoidable presence of magnetic impuri- 
ties with the concentration as low as 0.1 — 1 ppm. On a 
theoretical side it was shown^ that magnetic impurities 
at these concentrations can explain the puzzling anoma- 
lous energy dependence of the relaxation rate observed 
in diffusive Au and Cu wires£*&. 

It is well known that a significant concentration of 
magnetic impurities strongly affects conventional super- 
conductors!. Further, it was shown recently that even a 
small concentration of magnetic impurities leads to ex- 
ponentially small tails in the density of states within the 
gap &2ii2iiiii2. Because such subgap states can trap and 
release the quasiparticles, they can become an important 
for a number of physical applications, e.g. for Josephson- 
junction qubits. For example, since the resulting subgap 
states are localized, the quasiparticles that are gener- 
ated during qubit manipulations or readout — , can get 
trapped in these localized subgap states. These trapped 
quasiparticles provide extra degrees of freedom that con- 
tribute to the dephasing and dissipation in subsequent 
qubit operations. Although the typical materials used 
for the fabrication of Josephson-j unction qubits [ e.g., 
Al — AI2O3 — Al), are not intentionally doped with mag- 
netic impurities, some concentration of them is unavoid- 
able, especially in small devices where they are due to 
the surface effects. Therefore, it is important to estimate 
the effects of diluted magnetic impurities on the density 
of states in these systems and the resulting dephasing. 



With this goal in mind we present here a detailed the- 
ory of the subgap states in a diffusive superconductors 
containing low concentrations of magnetic impurities. In 
particular, we study the contribution to the subgap den- 
sity of states due to spatial fluctuations of the concen- 
tration of magnetic impurities. Combining the results of 
our analysis with those recently obtained by Lamacraft 
and Simons 12 , who considered the subgap tails resulting 
from universal mesoscopic gap fluctuations, we provide a 
full picture of the subgap states, and establish in which 
regime one or the other type of the fluctuation domi- 
nates the physics. Finally, the results of our theoretical 
analysis are applied to estimate the subgap density of 
states generated by a concentration of w 1 ppm of mag- 
netic impurities in aluminum-based qubits. The result of 
our analysis is that the smallest superconducting islands 
used in modern experiments might contain a few sub- 
gap states while larger qubits can contain about 10 3 of 
them. Finally, we give the estimates the qubit decoher- 
ence resulting from trapped quasiparticles and conclude 
that they provide an important source of the decoherence 
in larger devices but negligible in the smallest ones. 



The rest of paper is organized as follows. In SecHTIwe 
present the qualitative arguments and the results of our 
analysis. The technical details of the density of states 
computation are presented in Sec. IIIII there we first dis- 
cuss the effect of fluctuations on the disorder averaged 
density of states in a large collection of superconducting 
grains [ Sec. IIII Aj . and then extend the analysis to the 
bulk systems (d = 1,2,3) in Sec. llTTBl In Sec. ITVl we 
present the details of our estimates of the subgap density 
of states in superconducting qubits and its effect on the 
decoherence. Finally, Sec. IVl gives the conclusions. 
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II. OVERVIEW AND MAIN RESULTS. 



The effects of impurities on conventional superconduc- 
tivity were extensively studied in the last four decades. 
It is well known that dilute non-magnetic impurities pre- 
serving the time reversal symmetry do not affect the su- 
perconducting gapi^ while a finite concentration of mag- 
netic impurities suppresses both the critical temperature 
T c and the gap in the local density of states, E g . The 
effect of weak magnetic impurities is controlled by the 
dimensionless parameter 



TcA 



(1) 



where r s is the spin-scattering time and A the bare super- 
conducting gap. At the critical concentration of magnetic 
impurities corresponding to £ = 1 the gap in the density 
of states closes but T c remains finite; at larger concentra- 
tions one the superconductor becomes gaplessi Subse- 
quent extensions of this theory have included the effect of 
strong magnetic impurities. While a single strong mag- 
netic impurity generates a bound intra-gap quasi-particle 
state and a local suppression of the BCS gapi&*i£*i£, a fi- 
nite concentration of them leads to the formation of an 
intra-gap impurity bandi£*i2*i^. As the concentration is 
further increased, the impurity band broadens, merging 
with the BCS continuum, and eventually resulting in a 
complete gap destruction. 

For small impurity concentration corresponding to £ < 
1, the Abrikosov-Gorkov theory 7 predicts a vanishing 
density of states for | E \< E g = A(l - ( 2 / 3 ) 3 / 2 . This 
conclusion, however, neglects the effects of the impurity 
distribution fluctuations^. One expects that such fluctu- 
ations smear the gap edge and lead to the exponential 
tails in the density of states extending down into the su- 
perconducting gap, similarly to the Lifshits tails^i in the 
density of states below the band edge of a disordered 
conductor 21 i 22 i 23 . 

As in the case of Lifshits tails one can visualize the the 
subgap states in a bulk sample as appearing locally in 
places where the impurity potential decreases the energy 
of the quasiparticle. Further, the main contribution to 
the states deep in the tail is dominated by the specific 
impurity potential that is called the optimal fluctuation. 
Historically, the first example of such optimal fluctua- 
tions in the context of superconductivity was found by 
Larkin and OvchnikoA»2i who considered the effects of 
the interaction constant inhomogeneity on the density of 
the subgap states. Recently, various mechanisms lead- 
ing to such sub-gap states in superconductors containing 
magnetic impurities have been considered, both in lightly 
magnetically doped superconductors^, in small super- 
conducting grainsSiiiiiSi, and in diffusive superconductors 
We will concentrate on the latter case: dirty super- 
conductors containing both magnetic and non-magnetic 
impurities, under the conditions C 1 and I <C £, where 



I is the mean free path relative to momentum scattering 
and £ is the coherence length. 

The subgap tails in this regime were recently consid- 
ered by Lamacraft and Simons^. Slightly above the gap 
the density of states predicted by Abrikosov and Gorkov 



1_ j E-E g { Q 
TrL d * 



A- 3 / 2 ~ nv L 




(2) 
(3) 



where i>o is the density of states per unit volume in the 
normal state, L is the sample size, and d the dimension- 
ality of the system. Formulating the problem in terms of 
a supersymmetric Sigma model, Lamacraft and Simons 
obtained a density of states close but below the gap 



— cx exp 



where dd ~ I and 



"' / E - E^ 3/2 
A„ 



A 



E g -E 



1/4 



(4) 



(5) 



is the linear size of the optimal fluctuation. 

This result has a clear physical interpretation. In 
zero dimensions, the action is a universal function of the 
rescaled energy (E g — E)/A g , where A g can be inter- 
preted as the effective level spacing in the right above 
the AG gap. This indicates that the gap fluctuations can 
be seen as Random Matrix-like fluctuations of the edge 
of a Wigner semicircle 2 ^. Indeed, the zero dimensional 
result was previously conjecturedii on the basis of the 
universality hypothesis of Random Matrix Theory. 

In a bulk system [d ^ 0], the subgap localized states 
giving the density of states Eq.Q originate from par- 
ticular configurations of normal and magnetic impurities 
favoring mesoscopic fluctuations of the gap edge. On the 
other hand, it is intuitively clear that another type of in- 
homogeneous fluctuations should also contribute to the 
sub-gap tails of the DOS: local fluctuations of the con- 
centration ni mp of magnetic impurities dictated by their 
Poissonian statistics . A local increase of ni mp , and there- 
fore of the spin-scattering rate above its average value 
implies a local suppression of the AG gap, and as a result 
the formation of sub-gap localized quasi-particle states 2 ^. 

Below, we will show that the mechanism generating 
subgap states described above is complementary to uni- 
versal mesoscopic gap fluctuations. This complementar- 
ity is particularly transparent in the peculiar zero dimen- 
sional limit, where the explicit form of the subgap tails 
depends on the way the statistics over disorder realiza- 
tions is acquired . Indeed, if gap fluctuations are studied 
in a single sample, e.g varying the boundary conditions 
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through gate voltages, the rate £ is constant, and the uni- 
versal result of Lamacraft and Simons 12 always applies. 
On the other hand, if one considers a large collection of 
the superconducting grains C is g°i n g to display sample 
to sample fluctuations. As shown below, depending on 
the system parameters, either the former or the latter 
effect is bigger. 

While in the zero dimensional case one has the possi- 
bility to single out mesoscopic gap fluctuations by con- 
sidering a single mesoscopic grain, in finite dimensional 
systems this is not possible. Since the analysis of Ref. 0] 
neglects large scale fluctuations of limiting the scatter- 
ing rate to be a constant independent on position, only 
a direct comparison of the resulting subgap tails can es- 
tablish which of the two mechanisms gives the dominant 
contribution to the asymptotic subgap density of states. 

In the following Sections, we shall study in detail the 
physics of the subgap states due to the fluctuations of the 
concentration of magnetic impurities, or, equivalently, of 
C- We describe the dirty superconductor by the quasiclas- 
sical Usadel equations^ 2 ^ in which the spin-scattering 
rate £ becomes a position dependent statistical quantity. 
The variable £ inherits its statistics from the Poissonian 
distribution of magnetic impurities. In this framework, 
we find the optimal fluctuation of ( by solving the as- 
sociated variational problem and the resulting instanton 
equations, and discuss the deep analogies between the 
problem at hand and the Lifshits tails in disordered con- 
ductors. We then calculate the expression of the asymp- 
totic subgap tails of the DOS, including gaussian prefac- 
tors, as a function of the distance from the gap edge 



Se 



(Eg ~ E) 



and of the average scattering rate impurity concentra- 
tion rii mp , and dimensionality [d = 0,1,2,3]. For every 
dimensionality, we compare our results with the subgap 
tails due to mesoscopic gap fluctuations and establish the 
parameter that controls which of the two types of fluctu- 
ations dominated the physics of subgap states. 

The results of our analysis may be summarized as fol- 
lows. While the typical size of the optimal fluctuations 
associated to a local increase of £ is given by Eq. (J5J , the 
subgap density of states originating from the ni mp fluc- 
tuations is 



— oc exp 



exp 



-ad 



£4/3 ^ I 
\ \ / TP Z? x ^ 



L 



6E n 



(6) 



where ad is a constant, calculated below, 8E g = 
A( 2 / 3 /y/N, and N = n imp L d . A more detailed formula 
for the DOS is given by Ea. l|51[l . The effective dimen- 
sionality of the system is determined by comparing the 
linear size of the sample to the typical size of an opti- 
mal fluctuation A, Eq.JSJ). In particular, for d = the 



parameter that determines whether the physics of tail 
states is dominated by mesoscopic fluctuations [Eq. 
or fluctuations of £ [Eq. (©] is the ratio 



8E a 



For /3 > 1 the asymptotics is dictated by Eq.J3J, while 
for (3 <C 1 the result of Eq.© applies. In higher dimen- 
sions, a detailed comparison is more involved. However 
it is interesting to notice that, apart from a numerical 
constant, the ratio of the actions relative to the two opti- 
mal fluctuations [see Eq. J3J|, and Eq. @] is independent 
of dimensionality. 



III. SUBGAP DOS ASSOCIATED TO 
FLUCTUATIONS OF £. 



In this Section, we present the theory of the subgap 
density of states associated to fluctuations of the spin- 
scattering rate Cj starting with the simple zero dimen- 
sional case, and subsequently extending it to the case of 
finite dimensional systems. 



A. Fluctuations of £ in small grains. 



We begin with the simpler zero dimensional case, the 
results in this subsection apply to the samples which are 
smaller than the size of the optimal fluctuation. 

As explained in the previous section, the zero dimen- 
sional case is peculiar because the form of the subgap 
DOS depends on the whether statistics is acquired mea- 
suring the DOS in a single sample, e.g. varying the 
boundary conditions, or in a large collection of small 
grains. In the first case, the subgap tails are always 
dominated by mesoscopic Random Matrix-like gap fluc- 
tuations. In the following, we will analyze the case of a 
large collection of small grains, where both mesoscopic 
gap fluctuations and fluctuations of £ contribute to the 
subgap tails. 

For every value of the dimensionless spin-scattering 
rate £, and slightly above the gap edge, the DOS in the 
Abrikosov-Gorkov approximation is given by 



J_ lE-E g (Q 
irL 3 ' 



A" 3 / 2 ~ nv L 




(7) 
(8) 



where L is the linear size of each grain. We can estimate 
the DOS tail resulting from sample to sample fluctua- 
tions of the concentration of impurities by promoting £ 
to be a statistical variable, with average C and variance 



4 



( 2 /N, where N is the average number of paramagnetic 
impurities per sample. The average DOS is 

4 



2 

C^73 V 3 



d(6() j E~E g (C + 6C) c -^ 



where t]q — £ 2 /N . Since 

E g (C + SQ _E g (Q 6( 
A ~ A C 1/3 ' 

introducing Se — [E g (Q — E)/A, we obtain 



+ 00 



d{S() v / 5(-( 1 / 3 5e e"^T 



C /3 Se 



25 



3/4 



Here 



s = - 



2 C 4/3 ^ oe J 2 ^ SE g 



(9) 



(10) 



is the action, assumed to be S » 1, and 

X2/3 

fl£ fl = A^ 



AT 

is the typical scale of gap fluctuations. The final result is 



1 



exp 



1 AT 



2 C 4 / 3 



Let us compare now these results to the asymptotic 
DOS resulting from mesoscopic fluctuations of the gap 
edgei£. These fluctuations are characterized by the ac- 
tion 



_4 ( E g -E 



3/2 



As consequence of the universality of this result, this ac- 
tion is of order 1 at energies below the gap corresponding 
to the effective level spacing right above the AG gap edge, 
i.e. 



E g - U: 



A„ 



At the same time, the action of Eq. <|l(J|l is of order one 
the distance from the gap becomes of the order of the 
typical gap fluctuations 

E g — LO f» SEg. 

It follows that the parameter that roughly determines 
which one of the two mechanisms dominates is 



For (3 ^> 1 the asymptotic tails are dominated by meso- 
scopic fluctuations and the result obtained in Ref. 
applies. In contrast, when /3 <C 1 fluctuations of the im- 
purity concentration dominate the physics, and Eq.Q 
describe the asymptotic tails. 

In more detail, if (3 -C 1, as the energy E is decreased 
from Eg towards the Fermi level, beyond E g ~E w A g the 
mesoscopic gap fluctuations become exponentially rare. 
However, for A g < E g — E < SE g the system is still in the 
range of the typical gaussian fluctuations of E g associated 
to fluctuations of (, which dominate the physics. As the 
energy is decreased further, both asymptotic result are 
applicable, and to determine which wins, one has to com- 
pare the actions. The direct comparison shows that in 
the range SE g <C E g — E <C 5E g //3 3 the density of states 
is dominated by the action associated to optimal fluctua- 
tions of £. At the crossover point E g — E~ SE g / f3 3 , both 
actions are of the same order S ~ Ssusy — I//? 6 3> 1 
but typically at this point the density of states is negligi- 
ble. For example, for — 0.3, at the crossover S ~ 10 3 , 
and exp [-S] is practically zero. Therefore, we conclude 
that the simple rule to determine which mechanism dom- 
inates is to compare the typical gap fluctuations to the 
effective level spacing: the largest wins. 



B. Bulk Optimal fluctuations. 

In this section, we study the effect of spatial fluctu- 
ations of the spin-scattering rate in relation to subgap 
localized states. In particular, we will consider a dirty su- 
perconductor containing both normal and magnetic im- 
purities, under the following conditions 



I < £ C < 1 and 



of" » 1, 



(12) 



where I is the mean free paths relative to non-magnetic 
scattering, £ is the coherence length, d is the dimensional- 
ity of the problem, and ni mp is the average concentration 
of magnetic impurities in the system. This set of condi- 
tions describe a disordered superconductor containing a 
relatively large number of weak magnetic impurities. 

The condition I <C £ implies that the problem can be 
studied in the framework of the semiclassical approxi- 
mation 28 . In particular, parametrizing both the semi- 
classical Green's function g(r, e) and anomalous Green's 
function f(r, e) in terms of a phase 9 by 



ff(r,e) = cos[0(r,e)], 
/(r,e) = isin[0(r,e)], 



one gets the Usadel equatio: 



,27.28 



V 2 0(r) + iesin[0(r)] - cos[0(r)] 
Csin[0(r)]cos[0(r)] = 0. 



(13) 
(14) 



(15) 



5E n 



(11) 



Here the unit length is £ = y/D/2A, and D is the diffu- 
sion constant, and e = Ej A. 
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In this section, we will study small deviations from the 
solutions of the uniform Usadel equation 

C C 
(ie - | cos(0)) sin((9) - (1 + - sm(0)) cos(0) = 0. (16) 

Let us derive a few well known properties of such solu- 
tions, corresponding to the Abrikosov-Gorkov mean field 
theory 7 . Setting 6 — —tt/2 — i arctanh(w), we have 



sin(0) = 
cos(0) = — i 



1 



u 



(17) 
(18) 



Therefore, the solution of Eq. (|16[) are all those u such 
that 

1 



e = u 1-C 



From this equation we obtain that the gap edge is at 
eo = (1 — £ 2 / 3 ) 3 / 2 . In particular, the density of states per 
unit volume 

v = ^ o Re[cos(0)], 

is uniform and equal to zero for | e |< eo- At the gap 
edge the parameter u is equal to Uq = (1 — C 2 / 3 ) 1 / 2 . 
Correspondingly, one has the solution 9 edge = —tt/2 — 
i arctanh(ito). 

The Usadel equation is obtained neglecting large scale 
fluctuations of the impurity concentration. This assump- 
tion limits the scattering to be homogeneous in the ap- 
propriate long wavelength limit. On the other hand, in 
a realistic system the distribution of impurities is poisso- 
nian. This implies that the number of impurities in any 
finite volume element fluctuates, with a variance equal to 
the average. 

Following this observation, we introduce an effective 
theory describing the fluctuations of the concentration 
of magnetic impurities on length scales of the order, or 
larger, than the coherence length. In terms of the Usadel 
equation Eq. Ijl5|l £ becomes a position dependent random 
variable. As a consequence of the Poissonian statistics of 
the impurities, we have 



C(r) = C + 5C(r), 
(SO = 0, 



(19) 
(20) 



(<5C(r)5C(r')> = 8(t-t>), (21) 

where £ is the average, uniform, dimensionless, spin- 
scattering rate, and d is the dimensionality. 

Let us look for solutions of Eq. (|15fl at energies slightly 
below the gap (e = eo — Se) in the form 

0(r) = Oedge - i<f>(r)- 

Expanding Ea. (|15|l in <p, Se, and <5£, in the limit ( <C 1 
one obtains 



v 2 + ^ 1/3 ^ 2 - 



i 



i 



^173 { Se ~^m s<: 



(22) 



Rescaling again the length in units A = £ (2/3) 1 / 4 , and 
defining V = (3/2 C 2/3 ) 1/2 </>, we recast Eq. in the 
simpler form 



where 



V 2 V + ^ 2 = Se-f(v), 



/(r) = 71 



(23) 



(1/3 • 



In particular, 

(/(r)/(r')) = ^(r-rO, 



n 



C 



4/3 



^imp £ \ ^ / 



d/4\ 



(24) 



Let us split if) as ip — —x + iy. Then we have the 
system 



V 2 x + x 2 -y 2 = 8e-f(r), 
--V 2 y + xy = 0. 



(25) 
(26) 



Interestingly, this set of equations is analogous to the 
equations obtained by Larkin and Ovchiniko\t24, in the 
context of the study of gap smearing in inhomogeneous 
superconductors. 

Notice that the analytical properties of the quasi- 
classical Green's functions impose some constraints on 
the solutions of these equations. Indeed, the DOS is 



v(r)/v = Re[cos(0)] 



3 C 4/3 



(27) 



where vq is the bare DOS per unit volume at the Fermi 
level. Therefore, we must have y(r) > 0. 

Our aim is to evaluate the average DOS (v(t))/vq = 

\/2/3 C 4 / 3 (y( r )) at a distance Se below the average gap. 
In particular, 

. / w f 2>[/(r)l y(r,Se \ f(r)) eT^ S dr ^ 2 
(»(')> = - ~ \ r • (28) 



jv[f(v)] e -^I^mr 



For energies below, but not too close, to the average 
gap edge, the leading contribution to the DOS comes 
from exponentially rare fluctuations of /(r), describing a 
local increase of the impurity concentration, and there- 
fore a local suppression of the gap below its average value. 
In particular, we have to find the optimal fluctuation of 
/(r), i.e. the configuration of /(r) associated to the dom- 
inant contribution to (i/(r)). 

In order to be characterized by a finite action, an 
optimal fluctuation has to be such that /(r) — * as 
r = ! r |— > oo. The asymptotic behavior of any 
solution of Eas. (|25l26|l with f(r) describing an optimal 
fluctuation is therefore 



y(r) — > as r — > +oo, 
x(r) — > Voe as r — > +oo. 



(29) 
(30) 
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Using Ea. l|26|) . the second condition implies 



v2 y rr 

► V oe as r 

2y 



-co. 



(31) 



Since the system is diffusive, and scattering is isotropic, 
it is natural to assume a spherically symmetric optimal 
fluctuation, its scale being A, i.e. f(r) = / ( r / A ). 
Then y(r) — y(r/X) which implies 



2;/ 



1/A 



2y 

Thus Eq.® implies A cx (Se)- 1 / 4 . 

The problem above has many analogies to the problem 
of localized states in the Lifshits tails of a disordered 



systei 



,21.22,23 



[see Appendix A]. The analogy becomes 



clear if we substitute V(r) = x- 
as 



i8e and rewrite Eq. i|26|) 



-V 



V(r) 



y(r) = -VSey(r), (32) 
V(r) — > as r — > +oo. 



As in the Lifshits tails problem, we are essentially look- 
ing for a potential well V{r) that admits a ground state 
at energy —y/Se. A ground state condition is important 
because the nodes in the "wave function" y(r) are not 
consistent with the analytical properties of the quasi- 
classical Green's functions. Despite the close analogy, 
it is important to notice two differences with the Lifshits 
tails problem: (i) our wave function y(r) is a component 
of the Usadel phase directly proportional to the DOS, 
and (ii) the statistical weight of the potential in the 
standard Lifshits tails problem is S cx J dr (V(r)) 2 while 
here it is 



S = — 
2i] 



2n 



dr (f(r)f 

dr (y 2 +V 2 V -V 2 -2V8e'V) 2 , (33) 



where we used Ea. (|25|) to express /(r) in terms of V(r) 
and y. 

With the exponential accuracy one can neglect the y 2 
term in Ea. (|33|l because the integration over the fluc- 
tuations near the saddle point characterized by Ea. lJ52*)) 
leads to a non-zero contribution to the density of states 
Eq. (|28J) [see Appendix B] . To find the parametric depen- 
dence of the action, we note that the Schrodinger equa- 
tion Ea. H32|) implies that the optimal fluctuation has a 
scale A cx (i5e) -1 / 4 and an amplitude Vo(r) cx — \/Se . 
Writing V — VSe v(r/\) we get 



S 



V 2 V 



V 2 - 2VSeV 



-d/4 



d(r/A) 



(V 2 v 



2v)' 



— const- 



L imp 



^4/3 



(«5e) 



2-d/4 



(34) 



where r\ is given by Eq.J2U). 

We now find the exact shape of the optimal fluctuation 
and the resulting numerical constant in Eq. I|34|l from the 
Euler-Langrange equations. We have to look for a the 
most probable fir) which when inserted in Eq. l|25l26|l . 
produces a solution y(r), positive everywhere, and such 
that S7 2 y/(2y) — > \J~8e at infinity. This general strategy 
can be applied to the classical problem of Lifshits tails, 
as explained in Appendix.A. 

Expressing f(r) in terms of x and y, we get the dimen- 
sionless action 



S = 



— { dr [V 2 x 

2y ' 



Se\ 



(35) 
(36) 



where d is the dimensionality. This action admits a trivial 
stationary point corresponding to zero action or f(r) = 
in Eqs I|25l2t)|l characterized by 



d-l 



Se = 0, 



(37) 



where the prime indicates the derivative with respect to 
r. It is possible to show that this equation admits in- 
stanton solutions. However, these solutions have cither 
y{r) — at every point in space, or y(r) < some- 
where. The analytic properties of quasi-classical Green's 
functions indicate that the latter solutions are not al- 
lowed. On the other hand, the solutions with y(r) =0 do 
not contribute to the DOS in the present approximation 
scheme. However, these solutions become meaningful in 
the Sigma model approach to this problem; there they 
describe the saddle point corresponding to mesoscopic 
fluctuations generating subgap statesi2i22i2i. 

Let us now find the equations associated with the non- 
trivial saddle points. In spherical coordinates, the Eulcr- 
Lagrange equation corresponding to the action, Ea. (|35|l . 
is 

d 2 ^- 1 A) -(d-l)d r {r d - 2 A) 
-2ir d -U=0, (38) 

A[x, x', x", r] = x" + d —^x' - x 2 + Se. (39) 
r 

The problem of finding the instanton solutions of this 
equation describing the optimal fluctuations is consider- 
ably simplified by the observation that, in all dimensions 
[d = 1, 2, 3], the solutions of 



2x' 



(40) 



are also solutions of Ea. 1)38(1 . Therefore, the nontrivial 
saddle points^ are described by the system 



x » + x > _ ^ + Se = 0, 

r 



d-l 



y = xy. 



(41) 
(42) 
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It is interesting to notice that the equation above for the 
variable x(r) coincides with the equation describing the 
trivial saddle point, Ea. 1)37(1 . in d— 2 dimensions. In three 
dimensions we get 



(43) 
(44) 



x = VSe (1 - 3 (scch [r/A]) 2 ) , 
yo oc — — tanh [r/A] (sech [r/A]) 2 , 



A = 



-1/4 



(45) 



The optimal fluctuation of the dimensionless spin-flip 
scattering rate can be then evaluated using 



fo(r) = Se 



2yo 



V 2 2/o 

2yo 



Thus, in three dimensions, the action at the nontrivial 
saddle point is 



(46) 



4tt 


(C 


384 


2t) 




5 


63 


(4/3 


6e 5/4 



In one and two dimensions the solution of Eq. (|41I42|1 is 

not so straightforward. In order to estimate the action of 
the optimal fluctuation, we had to minimize the action S 
in terms of a variational ansatz y a ,b{r) depending on two 
free parameters. The asymptotic form of any variational 
ansatz y a ,b{f) must be 



2r/A 




(47) 
(48) 



in order to obtain to a finite action. The result of this 
calculation is 



S 

s 



14 
2.8 



'-imp 



(4/3 

^imp 
(4/3 



Se 3 / 2 (d = 2) 
Se 7 ' 4 (d = 1) 



(49) 
(50) 



Clearly (y(r)) oc exp[— S]. This equality must be sup- 
plemented by the calculation of the prefactor, i.e. gaus- 
sian fluctuations, using the standard technique due to 
Zittarz and Langer— . This calculation is reported in 
some detail in Appendix. B. The result is 



~ " C 4/3 



5e ad exp 



- "-tmpC r 2— ri/4 
' ad C 4/3 d€ 



(51) 



where the exponent is ay — l/8(d(10 — d) — 12). The 
values of the numerical coefficient of the action a,d, as 
well as the estimated value of the numerical constants b e i 



according to the calculation reported in Appendix. B, are 
summarized in the following table 



dimensions 


ad 


b d 


1 


2.8 


0.0007 


2 


14 


0.006 


3 


63 


0.2 



First of all, the result of Ea. (|51l) is asymptotic, and 
requires S > 1. This condition translates into Se 3> Seo, 
where 



Sed 



i c 



4/3 



(52) 



Since we assumed ^ <gC 1, ni mp ^ d ^S> 1 , one has Se -C 1. 

As expected, the final result given by (|5*TJl crosses over 
to the result of the zero dimensional calculation per- 
formed in Sec. IIII Al in the limit d — > 0. The relation 
between the d ^ and d = results becomes simpler 
when we rewrite Eq. l|51|l as 



oc exp 



(53) 



where A = ^/(^e) 1 / 4 [see Eq.©]. Though this way of 
writing the result makes the d — > limit transparent, it 
should be kept in mind that the action in finite dimen- 
sional systems does not depend on the sample size L but 
only on intensive quantities. 

Let us now compare the asymptotic tails resulting from 
fluctuations of (, Eq.JSTJ with those associated to meso- 
scopic gap fluctuations 1 ^. The latter give an asymptotic 
DOS 



V{) 

SsUSY 



oc exp [Ssusy] , 
/47ri/ A£ d 



ad 



(Se) 



3/2-d/A 



(54) 



V C 2/3 

where ad is a numerical constant. In one dimension it 
is given by a\ ~ 8 4 -\/24/5 ~ 4000. Introducing Ao = 
^/(Je) 1 / 4 , one may write 



Ss 



USY 



E n — E 



3/2 



with A g given by Eq. (JSJ . From this equation we conclude 
that the ratio between this action and the action of the 
optimal fluctuations of ( [Eq.JSSJ)] is independent on di- 
mensionality (apart from a numerical prefactor) . This is 
a direct consequence of the fact that the typical linear 
sizes of the optimal fluctuations associated to the two 
mechanisms are identical. 

In order to establish which fluctuation dominates the 
physics of the subgap tails, we first notice that Ssusy is 
of order of unity at Se ~ Se' d , where 



i c 



2/3 



a d 4tw A£ c 



Therefore, in a finite dimensional system, the parameter 
that determines which one of the two mechanisms domi- 
nates is 



Localized Subgap States 



a - s < 
Pd = j—- 



(55) 



As in the zero dimensional case, when fid 3> 1 the sub- 
gap tails are dominated by mesoscopic gap fluctuations, 
and the asymptotics of the subgap tails is described by 
Ea. (|54|l . In contrast, when <C 1 the physics is dom- 
inated by the fluctuations of the concentration of mag- 
netic impurities, and Eq.JSTJ applies. More precisely, as 
the energy is lowered from the gap edge, i.e. as Se in- 
creases, first the asymptotic result relative to mesoscopic 
gap fluctuations starts being applicable beyond Se ~ Se' d . 
However, for Se' d < Se < Sed the system is within the 
range of typical fluctuations of £, which dominate over 
the exponential tails associated to mesoscopic gap fluctu- 
ations. Increasing Se beyond Sed both asymptotic results 
are applicable and the two actions should be compared. 
In particular, 



S 



(KY 



Ssusy 
implying that for 



Se 1 ' 2 , 



Se d < Se < Se ( — 



the asymptotic tails are dominated by fluctuations of £. 
At the crossover, i.e. Se ~ (Sed)^~ /{Se' d )^~ , both ac- 
tions are 



S ~ Ssusy 



»> 1. 



This implies that, at the crossover, the density of states 
is already negligible. For example, for d = 2 and at the 
crossover point, S ~ (l//3d) 3 . Taking, f3 d = 0.1, one 
obtains at the crossover S ps 10 3 , implying that exp(— S) 
is practically zero. 



IV. SHALLOW QUASIPARTICLE TRAPS IN 
SUPERCONDUCTING QUBITS. 



In the previous sections we found that a low concen- 
tration of the paramagnetic impurities lead to a small 
but finite density of the localized subgap quasiparticle 
states. In this section, we discuss the effect of these 
states on small superconducting devices, especially the 
ones proposed for quantum computation. Because the 
density of these localized states is always very low, they 
do not affect much the macroscopic properties such as 
specific heat or Josephson current of big junctions. How- 
ever, even a small number of quasiparticles trapped in 




rRAP sc 




sc 





FIG. 1: A pictorial representation of a Josephson junc- 
tion composed of two superconducting electrodes separated 
by an oxide layer, the basic building block of superconduct- 
ing Qubits. In these system, quasiparticles are typically ex- 
cited in the Qubit duty cycle. Delocalized quasiparticle diffuse 
through the system towards the metallic quasiparticle traps. 
In turn, quasiparticles trapped in localized subgap states re- 
main within the system. During subsequent Qubit opera- 
tions, these localized quasiparticles represent active degrees 
of freedom that may contribute to decoherence. Our esti- 
mates indicate that, for a junction with aluminum electrodes 
of size ~ 10 x 10 x 0.1 /im 3 containing w 1 pprn of magnetic 
impurities, the average number of localized subgap states is 
about 10 2 - 10 3 per electrode. The phonon-limited lifetime 
of a quasiparticle in such states is expected to be ~ 100 s at 
T = 10 mK. On the other hand, we estimate the presence 
of 10 2 — 10 s localized quasiparticles coupled to phase fluctu- 
ations in the junction to be limit the Qubit quality factor to 
10 4 — 10 5 , for an operational frequency of 10 GHz. 



these states might affect the long time phase coherence 
in the Cooper pair boxe o 33 ; 35 ' 36 or flux qubits 3 ^, so we 
shall focus on such devices here. 

The building blocks of these circuits are typically 
Al - AI2O3 - Al or Nb - A1 2 3 - Nb Josephson junc- 
tions. When quasiparticles are generated, say, during 
the readout process, most of them quickly recombine and 
disappear in the condensate. In order to eliminate the 
remaining quasiparticles, superconducting circuits incor- 
porate normal metal leads acting as quasiparticle traps 
13 i 35 . At low temperatures, these leads can be viewed 
as the sinks for all the quasiparticles that may diffuse to 
them. However, a quasiparticle localized within a super- 
conducting electrode cannot diffuse to the quasiparticle 
traps; instead it remains localized inside the supercon- 
ducting circuit and represents an active degree of free- 
dom that may contribute to decoherence in subsequent 
operations of the qubit [see Fig. ^ . 

We begin with the estimate of the number of local- 
ized states in a typical experimental setup. We fo- 
cus on aluminum based qubits. The typical parame- 
ters for aluminum are A = 200 /xeV, D = 50 cm 2 /s, 
v F = 2 x 10 s cm/s, and vq ps 1.5 x 10 22 cm~ 3 eV _1 . 
This implies that the coherence length is £ ~ 0.1 /im, 
while the mean free path is I ps 0.01 /im -C £. Therefore, 
the system is diffusive, and the first condition for the ap- 
plicability of the results obtained in previous sections is 
satisfied [see Eq. <|12[) ]. Further, because the thickness of 
the samples is about or less than the coherence length, 
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they are effectively two dimensional. 

The strength of the spin-flip scattering and the number 
of paramagnetic impurities is more difficult to estimate 
reliably. The problem is that for a typical sample of 
thickness I w O.lfim surface scattering is rather impor- 
tant and it is very likely that there it contains a signifi- 
cant spin-flip contribution. Indeed, the surface of the Al 
devices is covered by a thin layer of glassy Al 2 3 . A typ- 
ical glass has a concentration of utls ~ 10 7 — 10 8 /zm -3 
of two level systems which are usually attributed to the 
trapped electrons^ 7 .. This translates into a surface den- 
sity of 10 3 — 10 4 ^m~ 2 of free spins in the boundary layer 
of AI2O3 that interact with the electrons in the metal or 
into the effective density 10 4 — 10 5 ^m -3 per unit vol- 
ume of the device. Because this estimate counts only the 
spins associated with the two level systems in a glass, we 
expect that the actual number is somewhat bigger. On 
the other hand, assuming that the strength of the spin- 
spin interaction J »s 0.1 eV and its range r m 4 A we get 
that the volume density n s /~10 5 /z77i -3 leads to a spin- 
scattering time of t s wl0~ 8 s which is the lower bound 
given by the weak localization magnetoresistance data 3 ^. 
Thus, we believe that the realistic estimate for the effec- 
tive density of the paramagnetic impurities in the wires 
of the 0.1/itm thickness is about n s /~10 5 /im~ 3 . 

For this value of the the paramater (3 [see Eq. iJSHJ] 
is very small: fa ~ 1CP 5 <C 1. Therefore, the asymp- 
totics described by Eq.JSlJ applies. Using now Eq.JSlJ 
for the density of states we estimate the total number 
of the localized states per unit area in these conditions: 
p ~ 1 — 10 /im~ 3 . This means that for the large junc- 
tions of the area A » 10 x 10 /1m 2 , one gets about 
N su b s» 10 2 — 10 3 localized states on each superconducting 
electrode. The length scale associated to these states is 
about 1/Ltm. With the parameters above, we expect these 
localized states to be very shallow, typically 0.01 \ieV 
below the gap edge. This energy scale is two orders of 
magnitude smaller than the typical temperature in qubit 
experiments, T ~ 10 mK ~ 1 peV/kB- This implies 
that if a non-equilibrium quasiparticle generated in the 
duty cycle of the qubit gets trapped into one of these 
localized states, it will eventually absorb a phonon, get 
excited above the mobility edge, and diffuse out of the 
sample. However, at such low temperatures, the time 
scale r tr associated to these phonon-assisted detrapping 
processes is very long, r tr ~ 100s 40 i 41 . At the same 
time, in a superconducting qubit composed of large junc- 
tions of area 10 x 10 \irr? and operating at a frequency 
/ = 10 GHz the qubit relaxation rate associated to the 
coupling between the phase degrees of freedom and a 
single localized quasiparticle is T qp w 10 3 Hz [see Ap- 
pendix C]. Therefore, the presence of 10 2 — 10 3 localized 
quasiparticles interacting with the phase degrees of free- 
dom across the junction is expected to limit the quality 
factor to be Q = vo/T qp as 10 4 — 10 5 . It is impor- 
tant to mention that for smaller systems, e.g. of area 
A w 1 x 1 /im 2 , the probability of having a single local- 
ized subgap state of area considerably smaller than A is 



negligible. The subgap tails are effectively zero dimen- 
sional and the tail states extend to the whole system. So, 
an electron trapped in such state can always escape to 
the lead through the Josephson contact. 



V. CONCLUSIONS. 



We have shown that in a diffusive superconductor 
containing magnetic impurities at small concentrations 
[C < 1], the subgap tails of the DOS have two distinct 
contributions, one due to the universal mesoscopic gap 
fluctuations^, the second resulting from inhomogeneous 
fluctuations of the spin scattering time £■ 

We calculated the contribution to the density of states 
due to the fluctuations of ( [see Eq.© for d = and 
Ea.Q51|) for d > 0] and established, through a direct com- 
parison of the results, the parameter that controls which 
of the two mechanisms dominates the physics of subgap 
states [see Ea.l|ll fl -l|55 )l ]. The two mechanism are related 
by dimensional reduction, both at the level of instan- 
ton equations and in the dependence of the actions on 
5e. The deep reason behind the occurrence of dimen- 
sional reduction in this context is presently unclear, and 
requires further study. 

On the theoretical side, we believe that the mechanism 
generating subgap localized states studied in this work is 
not limited to the dirty superconductors containing para- 
magnetic impurities. For example, our study is clearly 
related to the study of inhomogeneous superconductors 
performed in a seminal paper by Larkin and Ovchinikov, 
Ref. |24| . In particular, in three dimensions they obtained 
result similar to ^5ip. It is interesting to notice that a 
recent solution of the same problem by the replica Sigma 
model approach^ gave the same asymptotic behavior as 
the mesoscopic gap fluctuations [Eq.Q]. The source of 
this discrepancy was previously attributed to the impos- 
sibility to construct in the present context a Lifshits ar- 
gument 29 that was used in the estimate of the action 
of the optimal fluctuation in Ref. pij . In contrast, we 
have argued that both results correctly describe two dif- 
ferent physical effects that work in parallel: mesoscopic 
Random-Matrix like gap fluctuations 2 ^ and long wave- 
length fluctuations of the coarse grained gap^i. In anal- 
ogy with the present study, a direct comparison of the 
two resulting actions determines which mechanism dom- 
inates the asymptotics of the density of states. 

Finally, we applied the theoretical results summa- 
rized above to estimate the effects of a weak concentra- 
tion of magnetic impurities [~ 1 ppm, corresponding to 
£ rs 10~ 4 ] in Aluminum-based Qubits and Josephson- 
junction arrays [see Fig.JIJ]. Our estimates indicate that 
the presence of a small concentration of magnetic impu- 
rities is mostly relevant for large qubits, designed around 
the 10 fim scale. Under the conditions above, Cooper 
pair boxes of the size 10 x 10 x 0.1 /zm 3 are expected to 
have an average of w 10 2 — 10 3 two dimensional subgap 
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localized states per electrode. At a base temperature of 
10 mK, the lifetime of a non-equilibrium quasiparticle 
localized in one of such states is expected to be as long 
as 100s. We estimate that this would limit that quality 
factor of the qubit operating at the frequency / = lOGhz 
by Q fa 10 4 - 10 5 . 
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APPENDIX A: ANALOGIES TO THE PROBLEM 
OF BAND TAILS IN DISORDERED 
CONDUCTORS. 

In the following, we illustrate how the variational rea- 
soning employed in the main text works in the context 
of the problem of Lifshits tails, and how the exact so- 
lution in one dimension for the leading exponential de- 
pendence of the DOS tails relative to this problem is 
recovered2ii22i2^. 

Consider the Schrodinger equation 



2 



V(r) 



y (r) = -E y(r), 



where V(r) is a white noise potential 

(V(r)V(r')) =jS{r-r'). 



(Al) 



(A2) 



Well below the band edge (E = 0) , the leading contribu- 
tion to the DOS associated to Ea. (|Al|) comes from rare 
local fluctuations of V{r) able to generate bound states. 
In order to find the leading exponential dependence of 
the average DOS as a function of energy, we have to find 
the most probable potential V(r | E) admitting a bound 
state at energy — E, i.e the optimal fluctuation. The av- 
erage DOS is then 



(P) « e s , 
S = ±J dr (V(r | E))\ 



(A3) 



Let us now look for a ground state at energy —E, i.e. 
assume that the wave function y(r) does not have nodes 
[and can therefore be chosen to be positive]. For every 
given y(r), we can use the Schrodinger equation to obtain 



V(r) 



2y 



E. 



(A4) 



This equation specifies for every positive, non- vanishing, 
and smooth y(r), the potential well V(r) that admits it 



as a ground state of energy —E. Notice that y(r) must 
tend to as | r |— + +oo. The same is true for V(r), 
otherwise its action would be infinite. Thus 



2y 



E as 



+oo. 



Consider now the functional 

I[V(r)} = f dr (V{r)) 2 . 



(A5) 



(A6) 



Let use Eq. (|A4(I to rewrite this functional as a functional 
of y(r). One obtains 



Mr)] 



dr 



y 



2/y 



- E 



(A7) 



The problem of finding the most probable V(r) admit- 
ting a ground state at energy — E is now reduce to the 
problem of finding the stationary points of I[y{r)\. Since 
the argument of the integral is squared, one might be 
tempted to find the most obvious stationary point by 
setting 



V 2 y 

2y 



E = 0. 



(A8) 



This procedure is obviously incorrect, since it corre- 
sponds to the search a bound state in a flat potential 
landscape V(r) = 0. In mathematical terms, Ea. (|A8|) 
does not admit nontrivial solutions satisfying the appro- 
priate boundary conditions introduced above. 

Therefore, let us go back to Ea. (|A7|) and find the non- 
trivial saddle points. This program can be definitely com- 
pleted in one dimension. First of all set 



Obviously, 



V - = f + .f 

y 



(A9) 



(A10) 



Let us now rewrite the functional in Eq. i|A7|l in terms of 
/(r). We obtain 



I[f(r)} = 2 



dr [f + f-2E 



(All) 



Using the Euler-Lagrange equation, one may show that 
the nontrivial saddle point is given by the solution of 



/ - 2/ 3 + 4 £ / = 



(A12) 



This equation describes the motion in an inverted double 
well potential with maxima at f± = ±y/2E. The corre- 
sponding instanton is a trajectory going from / + to /_. 
Integrating this equation one obtains 



f(r) 



/ 2^tanh(%/2~Er) 



(A13) 
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Correspondingly we find 



Using Ea. (|Bl|) - (|B2(l we can estimate 



yo (r) oc sech(V2Ex) 



(A14) 



Finally let us evaluate the action corresponding to the 
potential relative to yo(r). We have to calculate S = 
l/(2^/)I[y Q (r)]. Inserting the expression for y and eval- 
uating the integral one obtains 



S=l-{2Ef' 2 
37 



(A15) 



This is corresponds to the action obtained by Halperin 
in its exact solution of the Id problem, Ref . [2lj . 



APPENDIX B: GAUSSIAN FLUCTUATIONS. 

In this appendix we outline the calculation of the full 
expression of the subgap tails of the DOS, including gaus- 
sian fluctuations. 

For every function f(r), we select the center of the r' 
of the optimal fluctuation in such a way as to minimize 
the functional 



D(f\r')= / dr[f(r)-f (r-r')} 



(Bl) 



Of all the solutions, we select the r' closest to r, the 
point at which the DOS is evaluated. Subsequently, the 
function f(r) is expanded in an complete set 



f(r) =X)f»^( r )' 



(B2) 



where the first d + 1 functions of the expansion are se- 
lected according to 



<h(r) = a fo(r - r'), 
Mr) = W/ (r-r'). 

Here the constant a and b are given by 

-1 1/2 



dr(Mr)Y 



~ I .IrlVMr))- 



1-1/2 



(B3) 
(B4) 



(B5) 
(B6) 



Having in mind this expansion, we write 
(y(r)) = ( f dr"y(f(r))6(r" r')> 
= < / dr"y(f(r))5(VD(f \ r")) \ detVVD |),(B7) 



where a is a region of space containing at most a sin- 
gle optimal fluctuation, y(f(r)) the appropriate solution 
of Eq. l|25|l - i|2t)[) relative to f(r), and in the last line we 
operated a change of variables from r" to D. 



det WD 



3_ 

b 2d- 



(B8) 
(B9) 



At the same time, while to lowest order y(f(r)) oc yo(r — 
r'), the coefficient of proportionality is set by fluctuations 
of / around f . Using Ea. lj25|) . Ea. lp3|) . together with 
Ea. (|B2l) . one obtains 



y(f(r)) « Ayo(r) 
A = 



Co - a' 



y J dr<j} Q (r)[y {r)} 2 
Now using the results above, it is easy to evaluate 

(y) ~ dr"y(f(r))d(^)) 



(BIO) 
(Bll) 



1 

(27T77) 3 / 2 6 d 

J ,+00 



(27T77) 1 / 2 



^4 yaz 



dor 



where 



(27T77) 2 



Jdry (r) 



(far 
e 2 i 



/dr-0 o (r-)[yoW] 



2 & d ' 



(B12) 



(B13) 



The last object we have to evaluate is the integral de- 
scribing the contribution to the DOS due to gaussian 
fluctuations of the height of the optimal fluctuation, i.e. 



I 



dWCo - a -1 e *r 



= 77 3/4 r(3/2)exp[-5/2]D_ 3/4 (A/2^) 

/ 1 \3/4 

c r) 3 ^T(3/2) \±j e~ s (B14) 

where S is the action at the saddle point calculated be- 
fore, -D_ 3 /4 is a parabolic cylinder function, and the last 
line is the lowest order asymptotic expansion for S ^> 1 . 

Using now Eq.(BI ^ -(E ^ -(BI3 fr . and Ea.l|27|l one ob- 
tains the final result Ea. (|51|l . 



APPENDIX C: ESTIMATE OF THE QUBIT 
RELAXATION RATE 



In this Appendix, we estimate roughly the Qubit re- 
laxation rate due to the interaction between a single lo- 
calized quasiparticle and the phase across the Josephson 
junction. In the following, we will have in mind large 
Josephson junctions of area A ss 10 x 10 /jm [see Fig. (Q] 
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with electrodes of thickness l z ~ 0.1 fxm, operating at 
a linear frequency of / = 10 GHz, corresponding to a 
circular frequency of loq ~ 6 10 10 s _1 . 

Let us consider the phase tp across the junction evolv- 
ing in time according to tp(t) = cos(woi). Correspond- 
ingly, the voltage V across the junction will evolve as 
V(t) — h/2e ip. For typical junctions, having oxide lay- 
ers of thickness 10 A, the potential drop is located 
mostly in a thin layer of width d ~ 1 A within the super- 
conducting electrodes. The rate at which a quasiparticle 
localized on a region of linear size A ~ 1 fim is excited 
due its interaction with this time dependent potential, 
can be estimated as 

^-TEl /*Th(r) eV(r) ^(r) | 2 

8{e f - a - fiwo), (CI) 

where an d ar e the wave functions and energies 
of the initial and final state respectively. Using the fact 
that the interaction is mostly concentrated close to the 
surface S of the electrode, we can write 



where 



/ 



dnpi(r)ipf{r) | 2 S(e f 



(C2) 



fcwo). (C3) 



Since the support of the wave function of the initial state 
is concentrated on an area 4~AxAft:lxl /im 2 , 7 can 
be roughly estimated having in mind a three dimensional 
diffusive quantum dot of dimensions A X A X l z as 



4ir 2 v a A 2 



vq A" J s 



drdr' [G(r, r', e)]_ [G(r', r, e + ^ )]_(C4) 



where z^q is the DOS per unit volume, and 



[G]- 



~i{G a - G r ), 



(C5) 



G r ' a being the QD's retarded/advanced Green's func- 
tions. 

Performing a disorder average^ on Eq. ( IC4(I , one ar- 
rives at 



qp 



huj [ — 



1 £ 



(C6) 



where D is the diffusion constant. Finally, inserting this 
expression into Eq. l|C2|) . and using the typical parame- 
ters reported above [Xp ~ 1 A, D w 50cm 2 /s, « 
1.5 10 22 cm -3 eV _1 ], one arrives at T qp w 10 3 Hz. 
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